1 



The shear viscosity of holographic superfluids 

Makoto Natsuume *) and Masahiro Ohta **) 

KEK Theory Center, Institute of Particle and Nuclear Studies, 
High Energy Accelerator Research Organization (KEK), 
1-1 Oho, Tsukuba, Ibaraki, 305-0801, Japan 
and 

Department of Particle and Nuclear Physics, 
The Graduate University for Advanced Studies (SOKENDAI), 
1-1 Oho, Tsukuba, Ibaraki, 305-0801, Japan 

o 

£\j ■ We study the ratio of the shear viscosity to the entropy density for various holographic 

superfluids. For the s-wave case, the ratio has the universal value l/(47r) as in various 
holographic models. For the p-wave case, there are two shear viscosity coefficients because 
of the anisotropic boundary spacetime, and one coefficient has the universal value. For the 
f I ' (p + ip)-wave case, the existing technique is not applicable since there is no tensor mode of 

metric perturbations which decouples from Yang-Mills perturbations. Our results indicate 
that the shear viscosity does not show a singular behavior at the critical point for holographic 
superfluids. 

Subject Index: 121 

<D , §1. Introduction 

In the studies of holographic superconductors/superfluids^^ (See, e.g., Refs.[8|)- 
[TO]) for reviews), one often uses numerical computations or some approximations. 
' This is because the holographic superfluids are Einstein-matter systems and it is in 

general difficult to solve such systems. One approximation often employed is the 
"probe approximation," where the backreaction of matter fields onto the geometry 
can be ignored. While the approximation is enough to see the phase transition and 
to compute properties such as the conductivity, gravitational properties of these 
systems, in particular analytic results are largely intact. It is desirable to obtain 
gravitational properties of these systems analytically. We investigate this issue in 
this paper. We study n/s, the ratio of the shear viscosity to the entropy density for 
^ \ holographic superfluids. 

$_i ■ Studying n/s is interesting from another point of view. According to the AdS /CFT 

duality, the ratio is universal, i.e., 

s An 

in the large- N limit. This holds for all known examples which have been studied 
(See, e.g., Ref . [TT]) and references therein). Even though there exists several argu- 
ments to support the universality,™"^ it is still unclear why the universality holds 
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microscopically and how generic the universality is. The holographic superfluids 
provide yet another example of the universality. 

The holographic superfluids exhibit a second-order phase transition. At high 
temperatures (normal phase), the bulk geometry is given by the standard Reissner- 
Nordstrom-AdS black hole. The shear viscosity for the Reissner-Nordstrom-AdS 
black hole has been computed in Refs. 16)- 19), and the universality has been shown 
already. Thus, we focus on the low temperature phase (superfluid phase). 

Technically, the universality largely depends on the following two properties of 
the bulk theory: 

1. One can use the Kubo formula to compute r\ and carry out the tensor mode 
computation of gravitational perturbations. There are no other fields which 
transform as a tensor even if matter fields are present. 

2. The entropy density is given by the Bekenstein-Hawking formula as long as the 
gravitational action takes the Einstein-Hilbert form. 

In this paper, we consider three class of holographic superfluids, the s-wave, p- 
wave, and (p+ip)-w&ve holographic superfluids in (d+l)-dimensional bulk spacetime. 
Our results are summarized as follows : 

(i) The s-wave holographic superfluids are described by Einstein-Maxwell-complex 
scalar systems In this case, the universality holds with a modification of 
the technique in Ref. 14). 

(ii) The p-wave holographic superfluids are described by Einstein- Yang-Mills sys- 
tems.® In this case, the Yang-Mills field breaks the SO(d — 1) rotational in- 
variance on the boundary theory, which has two implicat ions . First, the hydro- 
dynamic limit is not described by a single shear viscositv l**)| Second, for d = 3, 
one does not have a tensor mode which decouples from the Yang-Mills field. 
(Namely, item 1 of the above list fails.) As a result, the existing technique is 
not applicable. However, for d > 4, one has the SO{d — 2) invariance. In this 
case, a tensor mode exists, and the universality holds for the shear viscosity 
associated with the tensor mode. 

(iii) The (p + ip)-w&ve holographic superfluid is described by the same system as the 
p-wave holographic superfluid (with d = 3), but the symmetry breaking pattern 
is different.® Although the metric keeps the SO(2) invariance, the Yang-Mills 
field breaks the 50(2) invariance. As a result, there does not exist the tensor 
mode which decouples from Yang- Mills perturbations and the existing technique 
is not applicable. 

Our results indicate that the shear viscosity has no singular behavior across the 
phase transition for holographic superfluids (See Sec. 14. 2\) . 

The plan of this paper is as follows. In Sec. [21 we consider r//s for the s-wave 
holographic superfluids. In Sec. El we consider anisotropic holographic superfluids, 
the p-w&ve and {p + ip)-wave holographic superfluids. For the (p + ip)-wave case, 
we identify the Yang-Mills perturbations which couple to the would-be tensor mode 



*' Sec Rcf. 20) for the shear viscosity of another holographic superfluid. 
**' In the context of the AdS/CFT duality, anisotropic shear viscosities have been computed for 
the noncommutative Af — 4 plasmaPS 



The Shear Viscosity of Holographic Superfluids 



3 



of metric perturbations. In Sec. [U we discuss the implications of our results. We 
discuss the shear viscosity itself of holographic superfluids. For a p-wave superfluid, 
the entropy density s has been obtained, so we can discuss one of shear viscosity 
coefficients itself. The shear viscosity coefficient in the superfluid state is lower than 
the one in the (unstable) normal state. 



§2. The s-wave superfluids 



2.1. Background 

The s-wave holographic superfluids are described by Einstein-Maxwell-complex 
scalar system: 

0»i Xy /=gI R -2A- -K x (|^| 2 ) f mn f mn 



16TrG d+1 J "{ 4 

- K 2 (\&\ 2 ) \V M & - iqA M &\ 2 - V (|^| 2 ) j 

(2-1) 

with the ansatz 

d-l 

ds 2 d+l = -g tt (r)dt 2 + g xx (r) ^ dxj + g rr (r)dr 2 , (2-2) 

i=i 

A = A t {r)dt , (2-3) 
<F = <F(r) . (2-4) 

Here, capital Latin indices M, N, . . . run through bulk spacetime coordinates (t, Xi,r), 
where (t, Xj) are the boundary coordinates and r is the AdS radial coordinate. Greek 
indices fj,,u,... run though only the boundary coordinates. K\, K 2 and V are arbi- 
trary real functions of the scalar field. This action includes not only the conventional 
s-wave holographic superfluids^'^ but also generalized modelsP^I "O "VVe impose the 
regularity condition on the metric at the horizon r = r^: 

9tt ->■ Cf(r - r h ) , g xx ^c x , g rr -)• c r (r - r ft ) _1 . (2-5) 

These conditions fix the Hawking temperature and the entropy density of the bulk 
geometry: 

1 fc7 c (d ~ 1)/2 

t -tM' ; ^ 

The model exhibits a second-order phase transition. At high temperatures, the scalar 
field vanishes and one obtains the standard Reissner-Nordstrom-AdS black hole. 
But at low temperatures, the Reissner-Nordstrom-AdS black hole becomes unstable 
and is replaced by a charged black hole with a scalar "hair." 

This system is supposed to be dual to some kind of superconductors/superfluids. 
In fact, the low temperature phase shows the expected behavior for superconduc- 
tors/superfluids. As superconductors, one can see the divergence of the DC conduc- 
tivity, an energy gap proportional to the size of the condensate, and the holographic 
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London equationP^^'®'® As superfluids, one can see the existence of the second 
and fourth sounds J23>I2HJ Also, vortex solutions have been constructed.^"® 

2.2. 77/s 

Since we are interested in the viscosity, the main object to study is the boundary 
energy- momentum tensor. According to the standard AdS/CFT dictionary,®"® 
the bulk gravitational perturbations act as the source for the boundary energy- 
momentum tensor. Thus, our task amounts to solve the bulk gravitational pertur- 
bations. 

Consider the fluctuations of the energy-momentum tensor which behaves as 
e~ tut . The fluctuations are decomposed by the little group SO(d — 1) acting on X{ 
(i = 1, • • ■ , d — 1). The fluctuations are decomposed as the tensor mode, the vector 
mode ( "shear mode" ) , and the scalar mode ( "sound mode" ) . 

One can use various methods to compute the shear viscosity. Among them, the 
most powerful one is the Kubo formula method, which uses the tensor mode: 

r, = - lim -ImG% 12 {uj,d) , (2-7) 

w— >0 LO 

where G^ 12 (uj,0) is the retarded Green function for the tensor mode T 12 at zero 
spatial momentum: 

G^ 212 (w,0) = -i 1°° d d xe iwt 9(t) /[T 12 (t,f),T 12 (0,0)]\ . (2-8) 

J —00 

To obtain the retarded Green function, we consider homogeneous gravitational 
perturbations which take the form 

9MN = 9mN + h MN , (2-9) 

where g$N 18 ^ rie background metric (|2-2[) , In the Lorentzian prescription of the 
AdS/CFT duality pHJ the retarded Green function (|2-8p can be calculated from the 
tensor mode h\2- We expand the action in terms of (j)(t, r) := /i 1 2(t, r) up to quadratic 
order and use the Fourier transformation 

Ht,r) = j t^tj e~ i0Jt+i *' s f k (r)<f> (k). (2-10) 

The retarded Green function is obtained as follows: 

1. Solve the classical equation of motion for the field fk(r) with the ingoing- wave 
boundary condition at the horizon and fk(f) — >• 1 at the boundary. 

2. Substitute the classical solution into the action and represent the action in 
terms of the boundary value 0o- Only surface terms remain, and drop the 
contribution from the horizon. 

3. The retarded Green's function is given by the kernel of the on-shell action: 

Son-shell = ~\ I -^M-k)G R (k)Mk) (2-11) 
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where the on-shell action is defined as Son-shell = (S + Sqh + S^t.) Ion-shell- Sgh 
is the Gibbons-Hawking term to provide a correct variational problem for the 
background geometry. S c .t. is the counterterm to renormalize divergences in 
the classical action. 

Thus, the problem is to solve the equation of motion for the field under the 
appropriate boundary conditions. 

From Eq. (|2-ip . the action which is quadratic in (j) is 



1 



16vrG d+1 



1 

2 

8 r 



^ yxx 



(2-12) 



with the help of background equation of motions (See Appendix IA.2p . Because of 
the little group SO(d — 1) acting on Xi, the tensor mode of the metric perturbations 
decouples from the rest of perturbations: the other modes of the metric perturbations 
Hmn, the gauge field perturbations 8 Am and the scalar field perturbation 5^. Thus, 
they can be set to zero consistently. Since the background geometry must satisfy 
the stationary condition, we add the Gibbons-Hawking term 

Sgh = I d d x^2K , (2-13) 



16vrG d+ i 

where 7^ is the boundary induced metric, tim is the normal vector to the boundary 
and K = j^V^riv is the trace of the extrinsic curvature of the boundary. This 
provides surface terms 

1 1 1 1 ■ ^W<?"W 1 



GH 



d d x 



-.0, 



>9t 



(2-14) 



Therefore, the bare action is 
1 



(2) 



(S s + Sgh) 



l6iTG d+ i 
1 



d d+l : 



7(0) 



1 



(V 



+ 



d d x 



9xx 



2 

y/-l 



(0) 



2g, 



>9r 




(2-15) 

The action diverges as r — > 00, so the counterterms at the boundary must be added. 
We need only the gravitational counterterm in order to evaluate the retarded Green's 
function for the energy- momentum tensor. According to the holographic renormal- 
ization procedure, the counterterm is 



S 



1 



c.t. 



16vrG, 



d+l 



+ 



d d Xy 



"7 



2(d- 1) 
Z 



+ 



d-2 



(d-4)(d-2) 



d 



4(d-l) 



(2-16) 
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where L is the AdS radius and [7] is the Ricci tensor made from the induced 
metric jau- These terms largely depend on the spacetime dimensions^- However, in 
order to evaluate the shear viscosity, we need boundary terms only up to first order 
in uj: 0(uj 2 ) terms in the action do not contribute to the Kubo formula because of 
the to — > limit. So, only the first term in Eq. (|2- 16|) is important and it becomes 

{2) Sc.t. = -r^r- [ d d xj^)(d-l)4> 2 , (2-17) 
IbirGd+i Jr^oo v 

for the tensor mode perturbation. This term removes the divergences from the second 
term of Eq. (|2-15p . As a result, the renormalized action is 



ItotG^Xs. + SGR + Si 



c.t. 



—M~ k ) [ - l^-^-f-k (r)d r f k (r) | 



-Jrr 



(27r) druv '\ 2 g 

+ (terms which are proportional to the EOM) 
+ (contact terms) , 



(2-18) 



Here, we neglected the second derivative respect to t because it provide only O (a; 2 ) 
terms, "(contact terms)" provide contact terms in the Green function and have 
the form f-kfk- They will not affect the shear viscosity since they do not give an 
imaginary part of retarded Gr een' s function. The terms which give the imaginary 



part take the form like f-kd r fk 1^1- We will see this at the end of this section. With 
(12401) 

In order to find the on-shell action, we need to solve the equation of motion for 
fk(r): 



(g rr V~- 



9 



(0) 



ft + ^W/fc + v _ - / f h = , (2-19) 

9tt g rr y /_ g (0) 

where the long wavelength limit k — > is taken since C(|A;|) terms in the action don't 
contribute to the Kubo formula. The equations of motion can be solved as follows. 
First, solve this equation of motion near the horizon and impose the ingoing-wave 
boundary condition. Second, find the solution over the whole region in the bulk up 
to first order in to. Finally, match these solutions. 

First, consider the near- horizon limit of Eq. (|2T9|) . With asymptotics of the 
metric (|2-5p 

f k (r) ~ (r/r h - l) ±iW = exp [ ± in In [r/r H - 1] ] . (2-20) 

where to := oj/AttT is the rescaled dimensionless frequency. The ingoing- wave solu- 
tion is given by fk(r) = exp [—in In (r/rh — !)]■ We expand this solution in terms 



*' One has to be careful when the number of the boundary spacetime dimensions d is an even 
number. See Ref. 39) for details. 

**' So, the second term of Eq. (|2-15|) and the counterterm (|2-17p do not affect the shear viscosity. 
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of ro In (r/rh — 1) near the horizon since we take the ro — > limit at the end of the 
analysis. So, 

/ Jfc (r)~l-i»ln[r/r fc -l] . (2-21) 

is the boundary condition as r — > r^. The overall factor will be determined by the 
boundary condition at r — > oo. 

Next, get the solution of Eq. (|2-19p for all r. In order to evaluate the Kubo 
formula, it is enough to obtain fk{ r ) up to first order in ro. Thus, expand fk{ r ) hi 
power of ro: 

fk{r) = / (0) (r) + to/«(r) + O(tt) 2 ) . (2-22) 
Inserting this into the equation of motion, /v°) and /W satisfy 

(V r v /= <7^y 

/«" + ^ l==-f {i) ' = , (2-23) 

where i runs i = 0, 1. Solutions are given by 

/« ( r ) = C« + C 2 (i) /°° dr'-jtd=L= , (2-24) 

where C- s are integration constants. From the boundary condition at r — > oo, 

Cj 0) = 1 , = . (2-25) 

The rest of constants are determined by the boundary condition at the horizon. 
Since the integrand in Eq. (|2-24p has a simple pole at the horizon, 

dr > ^r(r') ^ [°° d(r>/r h ) = 1 _ 

>r~r h V-9 {0) (r') V ctct 1 J r /r h ~i W/r h - 1) l^G d+1 sT J ' 

(2-26) 

Comparing this with the boundary condition (|2-21|) . one gets 

Cf ] = 0, roC^ = 4G d+1 ■ ius . (2-27) 
Therefore, the solution of Eq. (|2-19|) with the appropriate boundary conditions is 

f k (r) =(l + AG d+1 ■ ius dr' ^Q— + O (ro 2 )") , (2-28) 



which becomes 



f k (r) -> 1, d r f k (r) -> -AG d+1 ■ icos grr(r) , (2-29) 

V-5 rlUj ( r ') 

as r — t- oo. Thus, the terms f-kfk and f-kdrfk hi the action provide real and 
imaginary parts, respectively. So, the contact terms, which have the form f-kfk: do 
not contribute to the shear viscosity. 
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Now, we are ready to extract the Green function from the on-shell action. Sub- 
stituting the solution (|2-28|) into the on-shell action, one gets 



Son-she,, = / ^- d M~ k ) (-l-^f- k (r)dMr) ] Mk) 



irr 



1 f d d k ~ 



iujs 



2 J (2ir) d ^ v ' \ Air 



This leads to the retarded Green's function 



G^(u,0) = -^ + o((u/T) 2 ) , (2-31) 

from the prescription f|2- 1 lj) . Finally, the Kubo formula (|2-7p derives the shear 
viscosity, 

^ = _ Hm 1^12^ * (2 . 32) 
s-0 CO 47T 

Thus, the ratio of the shear viscosity to the entropy density is 



(2-30) 



r\ 1 

S 4-7T 



(2-33) 



Seym = ... , , 
167rG rf4 



Therefore, the universality of 77/s holds in this system irrespective of whether the 
complex scalar condenses or not. 

§3. Anisotropic superfiuids 

The p-wave or the (p + ip)-w&ve holographic superfiuids are described by the 
Einstein- Yang-Mills system: 

1— J d d+l x^-g{R-2A-\{F a MN f} , (3-1) 

where F^ N = d M A a N - d N A a M + g YM e abc A b M A c N is the field strength of SU(2) gauge 
fields, <7ym is the Yang-Mills gauge coupling and e abc is the totally antisymmetric 
tensor with e 123 = 1. The gauge field is written as a matrix- valued form: 

A = A a M r a dx M , (3-2) 

where r a = a a /(2i) using the Pauli matrices, so [T a ,T b ] = e abc T c . 

3.1. The p-wave superfiuids 

The p-wave case is described by the ansatz 

d-l 

ds 2 d+1 = -g tt (r)dt 2 + g xlXl (r)dxj + g X2X2 (r) ^ dx\ + g rr (r)dr 2 , (3-3) 

i=2 

A = &(r)r 3 dt + wfr^dxi . (3-4) 
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The function <£(r) gives the background static electric potential whereas the function 
w(r) represents the condensate. We impose the regularity condition at the horizon 
r = r h : 

9tt -> Ct(r ~ r h ) , g XlXl -> c Xl , g X2X2 -> c X2 , # rr -> c r (r - r h )~ l . (3-5) 
Then, the temperature and the entropy density are given by 



1 [Cf V C X1 C X% 

T =7-\-> s = 1 Jr ' 3-6 

4vr V c r 4G rf+ i 

respectively. 

As is clear from the metric (|3-3p . the boundary spacetime is anisotropic. In such 
a case, the shear viscosity is no longer given by a single coefficient rj. Rather, one is 
interested in 

r] ijkl = - lim -lmG% kl (oj, 0) , (3-7) 
G% M (co,0) = -i d d xe iujt 6{t) /[T ij (t,x),T H (0,0)}\ . (3-8) 

From the symmetric nature of the metric and the SO(d — 2) invariance acting on 
%2,~" ,%d—i} there are only two nontrivial independent coefficients, e.g., rj 1212 and 
^2323 _ -^y e w ^y[ examine these coefficients below. 

The shear viscosities of anisotropic fluids have been widely discussed in the con- 
text of liquid crystal!^ There are various parametrizations known in the literature. 
Among them, the most well-studied parametrization is the Miesowicz viscosity coef- 
ficients. The coefficients rj 1212 and rj 2323 are related to the Miesowicz coefficients.® 
However, various conventions are found in the literature for the Miesowicz coeffi- 
cients. To avoid the confusion, we keep using rfi kl . 

3.1.1. r/ 2323 (for d > 4) 

First, let us consider rj 2323 . The coefficient exists for d > 4, and the metric 
has the SO(d — 2) invariance, so the perturbation /123 transforms as a tensor mode. 
Then, the discussion is similar to the s-wave superfluid case. 

The action (|3- 1 j) with appropriate boundary terms reduces to (h 2 3 =: (f>(r, t) = 

S^e lkx h(r)Mk)) 

167rG d+1 (2 \S p + S GH + S c . t .) 



d d k ( 1 ,OI \ 

X-*0 -7~—S-hif)MM Mk) 



Jrr 



+ (terms which are proportional to the EOM) 
+ (contact terms) , 



(3-9) 



using the ansatz (|3-3|) and the equation of motion for the background geometry (See 
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Appendix I A. 3, ip . The equation of motion is given by 



9 



(0) 



f h + + V ' '/'J fk = • (3-10) 

9tt g rr y /_ g (0) 

This takes the same form as the s-wave case (|2-19p . so the solution under the appro- 
priate boundary conditions is given by 

Afr)= { 1 + AGd+1 -^r dr ' J^(r') +0 ( ( " /T)2 )) ' (341) 

and the retarded Green function has the same form as the s-wave one (|2-3ip . From 
the Kubo formula, the shear viscosity to the entropy density ratio is 

„2323 i 

= J- • (3-12) 

s 4ir 

The universality holds for this CcLSG ELS well. 
3.1.2. n 1212 

Next, let us consider rj 1212 . For d = 3, this is the only shear viscosity coefficient. 
The perturbation h\2 transforms as a vector under SO(d — 2). Thus, it couples to 
the vector mode of the Yang-Mills perturbations. As a result, the existing technique 
is not applicable. 

It should be straightforward to obtain the action for the relevant vector mode 
perturbations since they are standard vector perturbations for which one can rely on 
the symmetry. However, it does not seem straightforward to solve them analytically. 
Thus, we will not discuss those perturbations further in this case. This is in contrast 
to the (p-\-ip)-w&ve case in next subsection, where it is not very clear how the relevant 
modes are coupled. This is because the symmetry structure is more complicated for 
the (p + ip)-w&ve case. 

3.2. The (p + ip)-wave superfluids 

For completeness, let us consider the {p + ip)-w&ve holographic superfluid. The 
(p + zp)-wave case is described by the ansatz 

ds 2 = -g tt (r)dt 2 + g xx {r)(dx\ + dx\) + g„(r)dr 2 , (3-13) 
A = $(r)r 3 dt + w^i^dxi + r 2 dx 2 ) ■ (3-14) 

As in previous models, the regularity condition at the horizon r = implies 

9tt -> c t (r - r h ) , g xx -)■ c x , g u ->■ c t {r - r^ 1 , (3-15) 

which leads to the temperature and the entropy density as 

1 [Cf c r 

(3-16) 



47T V C r 4Ga 
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respectively. It is argued that the (p + ip)-wave background is unstable and it turns 
into the p-wave background.® But the analysis was carried out only in the probe 
limit and the full analysis including the backreaction has not been done. 

Unlike the p-w&ve case, the metric is isotropic. Thus, the shear viscosity is 
described by a single coefficient. The anisotropy in the (xi, X2)-plane is caused by 
the Yang-Mills field. The condensation breaks the 50(2) rotational symmetry in the 
(xi, X2)-plane as well as the U(l) gauge symmetry. But, as is clear from Eq. (|3-14|) . 
it preserves a diagonal £7(1) which is a combination of the two. Thus, there does not 
exist the tensor mode which decouples from Yang-Mills perturbations. As a result, 
the existing technique is not applicable. 

Since the whole system does not have the 50(2) symmetry, it is worthwhile to see 
explicitly how Yang-Mills perturbations couple to the "tensor mode" perturbations. 
In Appendix IA.31 we derived the interaction of the Yang-Mills perturbations and 
the metric perturbations [Eq. (|A-29[) ] . For the tensor mode metric perturbations 
h 2 2 = — /i 1 i(:= 4>d) an d /i 1 2 = h 2 \{:= o d)H the interaction reduces to 

)h ij QH MN f% IN = h i j (F-fy i . (3-17) 



2 



Here, we defined 



Imn 
(F-f) 



D M a a N - D N a a M , (3-18) 
V M S ab + gYMe acb A c M , (3-19) 
F aiN tf N . (3-20) 

As is clear from Eq. (|3T7|) . the tensor mode h l j couples to (F ■ f) l j, which 
transforms as a tensor under the diagonal £7(1) symmetry, <j> d and (j) d couple to 

(F ■ f) d := (F ■ f) 2 2 = —(F ■ f)\ , (F ■ f) od := (F ■ f)\ = (F ■ f) 2 1 , (3-21) 

respectively. They contain the following components of 5A a M : 

a d := 5A 1 1 = -5A 2 2 , a od := 5A 1 2 = 5A 2 1 , (the other modes) = . (3-22) 

These perturbations (p d, <j>d, QW and are all coupled. The explicit form of 
(F ■ f) is given by 

(F ■ f) od = g rr g xx (d r w)(d r a od ) + g u g xx gYM$w(D t a d ), (3-23) 
(F ■ f) d = g rr g xx (d r w)(d r a d ) - g u g xx gYM<t>w(D t a od ), (3-24) 

which include the covariant derivatives of Oj: 

Dta d = d t a od + gYM&a d , D t a d = d t a d - gYM&a od . (3-25) 

Therefore, these forms mix a od and a d - 

Let us summarize how the tensor mode metric perturbations couple with Yang- 
Mills perturbations: 



*' In the s-wave and p-wave cases, the diagonal perturbation <j> d and the ofT-diagonal pertur- 
bation (f>d are completely decoupled. So, we have set (j>d = 0. But this does not hold for the 
(p + ip)-wave case as we will see below. 
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» The tensor mode metric perturbations 4> d and <f>d couple to the tensor (F ■ 
f)od and (F ■ /)<£, respectively, where (F ■ f) is made from the Yang-Mills 
perturbations. So, the action for (frod/d no longer takes the minimally-coupled 
scalar form. 

» The Yang-Mills perturbations a d/d couple to cj) d and 0<j. As a result, 4> Q d 
couples to (j)d through Oj. 

The complete action in terms of these "tensor mode" fluctuations are given by 
{2 \S (p+ip) + S GH ) = Y^G^ f dAx{2) + £gaugc + dm) ; (3-26) 

2 r 1 1 



2 L v 7 v ' J 2 

+ (surface term) 



(3-27) 



(2) £ g auge = V -5 (0) 5 XX £ [-<T(d, ai ) 2 + ^ M ^ 2 « 2 + 2"(A«i) 2 ] , (3-28) 



i=l 

2 



(2) A nt = y/-gWJ2M F 'fii ' ^ 



i=l 

where we defined two-component vectors as fa = (4>od,4>d), &i = (tW,ct<i), (-^ " /)i = 
((F ■ f)x, (F ■ 7)2) and « runs isotropic components i = 1,2. ^ 2 Omt is the interaction 
term we have discussed in Eq. (|3-17p . The mass-like function M(r) is defined by 

M(r) 2 := g rr g xx {d r wf + g xx g 2 M {g xx w 2 - g tt <P 2 ) w 2 . (3-30) 

The action leads to coupled equations of motion for 4>odi ^di a od and ad- It is 
difficult to solve them analytically, and it does not seem straightforward to obtain 
the shear viscosity to the entropy ratio. 



§4. Implications of the results 



We study rj/s for s-wave, (p + ip)-wave, and p-wave holographic superfluids. 
The shear viscosity for the s-wave superfluids satisfies the universality (|1T|) . The 
p-wave superfluids are anisotropic, and there are two nontrivial independent shear 
viscosities, r/ 2323 and 77 1212 . We show that one of the coefficients r/ 2323 satisfies the 
universality. 

On the other hand, for another coefficient ry 1212 of the p-wave superfluids and for 
the shear viscosity of the (p + ip)-w&ve superfluid, the gravitational perturbations in 
question couple to the Yang-Mills perturbations even in Kubo-formula method, and 
the existing technique is not applicable. We extract the modes which couple to the 
gravitational perturbations. For the (p + ip)-wave case, we write down the perturbed 
action for those modes. 

It would be interesting to solve the equations of motion to obtain r/ / s for those 
cases. If it turns out that they also satisfy the universality, these shear viscosities 
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will give highly nontrivial tests for the universality. If it turns out that they do not, 
these shear viscosities will give interesting counterexamples against the universality. 
The results are interesting in either way. 

There is another technique to derive rj/s in the membrane paradigm context.^ 
According to the method, transport coefficients in the boundary field theory can be 
determined by (i) the flow equation for r-dependent transport coefficients, e.g., the 
shear viscosity rj(r) and by (ii) their values at the horizon. If the tensor mode metric 
perturbation is written as a free scalar, the flow equation becomes trivial in the 
hydrodynamic limit: d r r](r) = 0. In this case, (rj/s) (boundary = (v/ s ) I horizon = l/4vr. 
This method works for rj in the s-wave case and for rj 2323 in the p-wave case. But 
it does not work for r] 1212 in the p-wave case and for rj in the (p + ip)-wave case. 
This is because the interactions of the metric and Yang-Mills perturbations provide 
a non-trivial flow equation d r rj(r) ^ 0. 

We now discuss the shear viscosity itself of holographic superfluids below. 

4.1. Viscosity of superfluids 

The holographic superfluid shows a nonzero viscosity. To interpret the result, 
note the following points. 

First, a superfluid has a nonzero viscosity. For example, for He no viscous 
resistance is observed when it goes through a narrow pipe, but a viscous drag is 
observed when a test body is moved in the liquid. 

According to the two-fluid model, a superfluid consists of the superfluid compo- 
nent and the normal component. The normal component has a nonzero viscosity, so 
a superfluid has a nonzero viscosity as a whole. The normal component represents 
the effect of thermal fluctuation, and it always exists at finite temperatures. And the 
quasi-particle description is valid for the normal component. Since we do not sep- 
arate the normal and superfluid components, one cannot observe the zero viscosity 
for the superfluid component. 

Second, currently the boundary theory description is not clear for holographic 
superfluids, but the boundary theory presumably contains the fields which may not 
play an important role in the superfluid behavior. Among the other things, the 
boundary theory should include the SU(N) non-Abelian gauge field, which is unlikely 
to play an important role. The computation of r]/s includes the dissipation not only 
from the normal component but also from these fields. 

Obviously, while r]/s is the same in both phases, r\ itself can have different 
functional forms. This requires the knowledge of s, and it would be interesting to 
compute it. In the probe limit, both phases are described by the same bulk geometry 
since the backreaction of matter fields onto the geometry is ignored. Thus, one needs 
the fully backreacted metric to find a nontrivial behavior. In particular, it would be 
interesting to see if rj in the superfluid state is lower than rj in the (unstable) normal 
state. Again one needs a fully backreacted metric, but analysis near the critical point 
or a numerical computation would suffice for the purpose. 

In fact, the entropy density s has been obtained for a limited class of holographic 
superfluids. Especially, Ref. |4"2]) obtained s for the (4 + l)-dimensional p-wave super- 
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fluid in the grand canonical ensembleEE They use the parameter a := K^/g. In our 
notation, a oc 1/flYMi and a. — > corresponds to the probe limit. Once the backre- 
action is taken into account, the p-wave superfluid undergoes the second-order phase 
transition only when a < a c , where a c ~ 0.365, and it undergoes the first-order 
phase transition when a > a c . Namely, the phase transition becomes first-order 
when the backreaction becomes large. 

According to their computation, s in the superfluid state is lower than s in 
the unstable normal state below T < T c at fixed chemical potential fj, = A^. See 
Fig. 3(b) of Ref. I4"2j) . Since j^ 2323 satisfies the universality, r/ 2323 in the superfluid 
state is lower than the normal state one. This may be due to the zero viscosity of 
the superfluid component. Needless to say, this statement is only for one coefficient 
of shear viscosities of one p-wave superfluid. At this moment, it is not clear if the 
same holds in general. 

4.2. Implication to dynamic critical phenomena 

We found that the universality oirj/s holds both for high temperature phase and 
for low temperature phase. In the second-order phase transition, critical phenomena 
occur and one has singular behaviors in physical quantities. In the dynamic case, one 
has singular behaviors in various transport coefficients P- 4 ' But our results indicate 
that there is no divergence in the shear viscosity. (Since the entropy density is the 
first derivative of the free energy, it is continuous across the phase transition. Thus, 
the universality of rj/s implies that the shear viscosity is also continuous across the 
transition.) 

More precisely, in the dynamic critical phenomena, the relaxation time of the 
order parameter diverges, which is known as the critical slowing down. In fact, for 
s-wave holographic superfluids, the relaxation time of the order parameter diverges 
near the critical point F» In general, when a system has a conserved charge, the 
associated transport coefficient diverges as well. For example, for T^ u , one has a 
(mild) singularity in n. But our results indicate that this does not happen in the 
holographic superfluids. The fact that singular behavior does not occur in rj has 
been observed in the critical phenomena of R-charged black holes P2> 
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For a s-wave superfluid, s has been obtained in the microcanonical ensemble.^ 
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Appendix A 

Quadratic forms of perturbations for Einstein- Matter actions 



In this Appendix, we derive the quadratic forms of the tensor mode perturbation 
for the s-wave, p-w&ve and (p + ip)-wave holographic superfluids. First, we derive the 
quadratic form of the Einstein-Hilbert action. Then, we derive quadratic forms of the 
matter action for these three models. Since we will not consider scalar perturbations, 
we will focus on the metric perturbations and the gauge field perturbations. 



A.l. The quadratic form of the Einstein-Hilbert action 



(0) . 



Consider the general perturbation Hmn to the background metric g MN - 



9MN = 9 M \ + h-MN ■ 



Under the perturbation, 
and 



± + \h+\(\h 2 -\h MN h MN 



(A-l) 
(A-2) 



(2)fl =V M (h IJ V M hjj + h MN V N h - h^V^N - h IJ V!hj M ) 



- \{V N h IJ )V N h u + \{V N h IJ )V!h JN - \v N hV N h + WR MI h T N h MN . 

(A-3) 



Therefore, the quadratic form of Hmn in the Einstein-Hilbert action is 
(2 ^f = " \h MN h MN ) {0) R + (h MI hj N - \hh MN ^ ^R MN + V M J M 

- \(V M h IJ )(V M hu) + ^(V M ^ J )(V^ JM ) + ^(V//i)(V J /0 - l -(Vjh IJ )(Vjh) , 

(A-4) 

where 

jM = h U V M hij + h MN VNh _ h MN Vih I N _ h U Vlhj M _ 1 h + ^ Wjv ^ . 

(A-5) 

If we restrict the perturbation Hmn to the tensor mode perturbation h\2, this 
quadratic form reduces to (4> := h l 2) 



(2) 



(A-6) 



where is the quadratic form of the Ricci scalar with the tensor mode perturba- 
tion: 



+ d r 2 



7^ 



(/><9 r </> + 



flrr 2<7; 



12^2 i,2 



ft 2 



(A-7) 



2:2^2 
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Here, we take the (d + l)-dimensional p-wave metric (|3-3p to preserve generality. 
Except the free scalar part, all these terms will be removed in the end. 

A. 2. The quadratic form of the s-wave holographic superfluid action 
The s-wave holographic superfluid is described by Eq. (\'2-l\) : 

S ° = TF^ f d d+1 x(^g~R + £,_ matt er) ; (A-8) 

IbirGd+i J 

1 



•^s-matter • — V 9 



K X F MN F MN - K 2 (|#|) \D M V\ 2 - V (|#|) 



,(A-9) 



where we defined a covariant derivative as Dm := Vm — iqA-u- Under the general 
gravitational perturbation 

9MN = 9m N + flMN , (A- 10) 

one can easily find 

(2) £ S - m attcr = \sf^h M N [ftX™" + K 2 Y MNIJ + V P MNIJ ] hjj , (A-ll) 

where 



P 



MN1J 



I (g{0)MI g (0)NJ + g (0)MJ g (0)NI _ g (0)MN g (0)IJ^ ? ( A . 12 ) 



X MNIJ ._ If ab F ab P MNIJ + - F A M F AN g^ IJ 

-F A M F AJ g^ NI - l -F MI F NJ , (A-13) 
yMNIJ ._ \ Da ^ P MNIJ + ( D M^! D Nyy g {0)IJ 

-2(D M &)(D J V)*gW NI . (A-14) 

Note that some modes of the metric perturbations couple with the gauge field per- 
turbations 5 A] and the complex scalar perturbation 5& in general. However, we drop 
these perturbations since these decouple from the tensor mode metric perturbations. 
The equation of motion for the background field is 

Q(0)^(0)M7V _ (0) R MN^ + T (0)MN = q (a . 15) 

The background energy-momentum tensor T$ N is defined as 

rp(0)MN _ 1 <9£ S -mattcr 

yJ-gW dg M N 
= \K, {-\g^ MN F ab F ab + F^ a F n A 

+ K 2 (-~g(°) MN (D A V)(D A &)* + (D^)(D N ^)A - ± g (°) MN v . 

(A-16) 



*> Here, we defined the symmetric symbol as F iM A F N)A = \ (F M a F na + F N A F A 
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This equation of motion leads to a relation between the Ricci scalar and the matter 
fields 

{0)R = 'a7T~~T\ Kl f mnF MN + K 2 \D M V\* + ±±±V , (A-17) 
4(a — 1) d—1 

and an isotropic component leads to 

R\ = -^l(v- I^FmkF^ . (A-18) 

So far our discussion does not assume an explicit background nor perturbations. 
Now, we take the s-wave background ansatz (|2-2p - (|2-4p and the tensor mode h \ = 
(j)(t,r). The quadratic form (|A-lip reduces to 

{2) C M ter = (hc^MNF 1 *" + K 2 \D M *\ 2 + V^ 2 , (A49) 

where we set the gauge field perturbations au to zero since these decouple from <j). 
Using the trace of the equations of motion (|A-17p and an isotropic component of 
Eq. (fA~48|) . one gets 

(2) £s- m atter = yf^P) (htiR - (°)i? 2 2 ) 4> 2 ■ (A-20) 

Combining the Einstein-Hilbert term (|A-4p and the matter term (|A-20p . we obtain 
the quadratic form of the s-wave holographic superfluid action (|2-12p : 



- \V=9(Vm<P) 2 - d t (2^gV)g U (t>d t (j> 
+ d r \^gW (2g rr 4>dA+\—9 rr 4> 2 

\_ \t Q X £C 



(A-21) 



with g xlXl = g X2 x 2 = 9xx- The second term (the total derivative with respect to t) 
does not affect the correlator, so we ignored the term in Eq. (|2-12p . 

A. 3. The quadratic form of the Einstein- Yang- Mills action 
The Einstein- Yang-Mills action is 

Seym = t-^ / d d+1 x(^=g~R + £ E YM-matter) , (A-22) 

167rG d+ i J 



1 771a zpaMN a 



(A-23) 



^EYM- matter 

where A a M is SU (2) gauge field and the field strength is defined as 

F a MN = d M A a N - d N A a M + g YM e abc A b M A c N . (A-24) 
Under the metric and gauge field perturbations 

9MN = 9 ( $ N + h MN , A a M = A(% + a% , (A-25) 



18 



Makoto Natsuume and Masahiro Ohta 



one can find*- 1 1 
^£eym 

(2V 



-grav 



(2)£ _l (2)£ , (2) £ . 

'—grav i '—gauge t *-mt 



1 



(2)£ 



gauge 



-5 



(o) 



1 



5 (0) 

Jmnj 

MNIJ f a 



-jhMN (X 
1 



MNIJ 



+ 2AP 



MNIJ 



)hu 



4 



aMN , x „a 7 MN b 

+ 2 M a6 ^ 



where 



/mjv := 

MNIJ ._ 



MNIJ 



yMN 
^ab 



Dmcl n — Dnci m 

1 77a EiaAB tjMNIJ , 1 P a AT rpaAN (0)1 J 

1 b AB b F +-t A b g y 

_pa MpaAJ J0)NI _ \paMN palj 

a y 2 , 

2 Q (0)MI paNJ _ 1 (0)MN paU 

y 2 y > 

,_abc jpcMN 



(A-26) 

, (A-27) 
(A-28) 
(A-29) 

(A-30) 
(A-31) 

(A-32) 
(A-33) 



The backgr oun d satisfies the Einstein equation (|A-15p with the energy-momentum 



tensor given by 

p(0)MN _ 



I fA g (0)MN pa jF aIJ _ p a(M A p\a\N)A _ g (0)MN 2A ) _ (A . 34) 



This equation of motion leads to a relation between the Ricci scalar and the matter 
fields 

(0) R d - 3 



4(d-l 

and an isotropic component leads to 



J7 a TP 



aMN + ^±1 2A 



d-1 



1 



d-1 



2A 



1 



T? a TP' 
^MN^ 



■aMN 



+ 



—F MX2 F« M *A 



(A-35) 



(A-36) 



A.3.1. The p-wave holographic superfluid action (tensor mode) 

Here, we derive the effective action of the tensor mode metric perturbation for 
the [d + l)-dimensional p-wave system. If we set metric perturbation (|A-25p to 
the tensor mode h 2 ^ = <f>(t,r), all the other perturbations are decoupled, so these 
perturbations can be ignored consistently. The quadratic action (|A-26j) reduces to 



(2) £ 



EYM-mattcr 



i(0) 



b MN b 



aMN 



+ 2A 



- g (0) fh0) R _(0) R 2\ ^ 

(A-37) 



*^ The quadratic form of the Einstein- Yang-Mills action was obtained in Ref. I47[) in order to 
calculate the one-loop divergence, but they omitted surface terms. 

**' Here, the vertical bars indicate that we do not symmetrize over a: F a< - M a F' a ' N ' A = 

1 tpaM A paNA + paN ^paMA^ 
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using Eqs. (|A-35j) and (|A-36j) . Then, one obtains the quadratic form of the p-wave 
action for the tensor mode metric perturbation: 

{2) S P = Y^G^{ I d d+1 x\f^ ~ \^M^f ~ d t Uyf^ftfkit, 

IV ^ 9x2X2 

(A-38) 

A. 3. 2. The (p + ip)-w&ve holographic superfluid action 

Let us derive the effective action of the "tensor mode" metric perturbations 
for the 4-dimensional (p + ip)-w&ve system. In this case, we must turn on four 
perturbations 4> od = h l 2 = h 2 ±, <p d = h}\ = —h 2 2, a 0( i = a 1 2 = o?\ and a d = a 1 ! = 
— a 2 2- The perturbed action is obtained by substituting these perturbations into 
Eq. fjA^H]) : 



(2)c 

° (p+tp) 



16vrG 



l — J d'x^gW (( 2 )£ grav + ( 2 )£ gaugc + ( 2 )£ int ) ; (A-39) 



{2) £ gra v = y£g® E \~\ {-9 tt (d t ^ + g"-(d r ^f} - l -M{r?tf 



8=1 

2 



,(A-40) 



(2) £gau gc = V -9 {0) 9 XX E [-^(dr^f + 4m» 2 «» 2 + 9 tt {D tai f] , (A-41) 
/ 2 

(2) Ant = V-5 (0) E^( F -^ > (A-42) 
i=i 

where we defined two-component vectors <fii = (<t>od>4'd), a i = ( a od,a>d)i (F ■ f)t = 
((F • f)x, (F ■ 7)2) and i runs isotropic components i = 1, 2. Here, we have omitted 
the surface term in Eq. (|A-40|) . The mass-like function M(r) is defined by 

M(r) 2 := /V* (M 2 + 9 XX 9ym [sT** ~ 9^) w 2 , (A-43) 
and the explicit form of (F ■ f)i is 

(F • f) od = g rr g xx (d r w)(d r a od ) + g u g xx g YM <Pw(D t a d ), (A-44) 
(i 7 • /) d = <T<f x (d^)(d r a d ) - <?!f° B !mi&v,{Pta od ). (A-45) 

Note that the covariant derivatives of have forms 

D t a od = d t a od + gvM^a d , D t a d = 5 4 a d - g Y M&a od . (A-46) 

Therefore, these forms mix a ot f and a^. 
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